Model Overview

We estimate the age of a site by calculating the years since the last fire. We then fit a curve to model the recovery of vegetation (measured using NDVI) as a function of it’s age. An additional level models the parameters of the negative exponential curve as a function of environmental variables. This means that sites with similar environmental conditions should have similar recovery curves.

Input data

The model was last fit on 2022-09-02 07:29:03.

This version of the model was fit with 108689 pixels including data from 2000-10-31 to 2014-06-26.

Workflow

This repository was developed using the Targets framework.

── Attaching packages ─────────────────────────────────────── tidyverse 1.3.2 ── ✔ ggplot2 3.3.6 ✔ purrr 0.3.4 ✔ tibble 3.1.8 ✔ dplyr 1.0.9 ✔ tidyr 1.2.0 ✔ stringr 1.4.1 ✔ readr 2.1.2 ✔ forcats 0.5.2 ── Conflicts ────────────────────────────────────────── tidyverse_conflicts() ── ✖ dplyr::filter() masks stats::filter() ✖ dplyr::lag() masks stats::lag()

Attaching package: ‘arrow’

The following object is masked from ‘package:utils’:

timestamp

Attaching package: ‘plotly’

The following object is masked from ‘package:arrow’:

schema

The following object is masked from ‘package:ggplot2’:

last_plot

The following object is masked from ‘package:stats’:

filter

The following object is masked from ‘package:graphics’:

layout

Registered S3 method overwritten by ‘hoardr’: method from print.cache_info httr Loading required package: abind Loading required package: sf Linking to GEOS 3.8.0, GDAL 3.0.4, PROJ 6.3.1; sf_use_s2() is TRUE

Attaching package: ‘lubridate’

The following object is masked from ‘package:arrow’:

duration

The following objects are masked from ‘package:base’:

date, intersect, setdiff, union

CmdStan path set to: /home/rstudio/.cmdstanr/cmdstan-2.30.1

Results

Environmental Drivers

These parameters represent the relationship of the following environmental variables to the recovery trajectory.

Recovery Trajectories

The plot below illustrates some example recovery trajectories. It currently just shows the top 20 cells with the most observations.

Model Performance

Compare estimated vs observed values for all pixels. This is not true validation - these pixels were included in the model fitting.

Spatial Predictions

Maps of spatial parameters in the model.

Park Reports

See the links below for park-level reports on vegetation status.

Addo-Elephant_National_Park Agulhas_National_Park Bontebok_National_Park Garden_Route_National_Park Karoo_National_Park Namaqua_National_Park Richtersveld_National_Park Table_Mountain_National_Park Tankwa-Karoo_National_Park West_Coast_National_Park